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ABSTRACT 

We report radial velocity measurements of the G-type subgiants 24 Sextanis (=HD 90043) and 
HD 200964. Both are massive, evolved stars that exhibit periodic variations due to the presence of a 
pair of Jovian planets. Photometric monitoring with the T12 0.80 m APT at Fairborn Observatory 
demonstrates both stars to be constant in brightness to < 0.002 mag, thus strengthening the planetary 
interpretation of the radial velocity variations. 24 Sexb,c have orbital periods of 452.8 days and 
883.0 days, corresponding to semimajor axes 1.333 AU and 2.08 AU, and minimum masses 1.99 Mj up 
and 0.86 Mj up , assuming a stellar mass M* = 1.54 M . HD 200964b, c have orbital periods of 
613.8 days and 825.0 days, corresponding to semimajor axes 1.601 AU and 1.95 AU, and minimum 
masses 1.99 Afj up and 0.90 M j up , assuming M* = 1.44 Mq. We also carry out dynamical simulations to 
properly account for gravitational interactions between the planets. Most, if not all, of the dynamically 
stable solutions include crossing orbits, suggesting that each system is locked in a mean motion 
resonance that prevents close encounters and provides long-term stability. The planets in the 24 Sex 
system likely have a period ratio near 2:1, while the HD 200964 system is even more tightly packed 
with a period ratio close to 4:3. However, we caution that further radial velocity observations and 
more detailed dynamical modelling will be required to provide definitive and unique orbital solutions 
for both cases, and to determine whether the two systems are truly resonant. 

Subject headings: techniques: radial velocities — planetary systems: formation — stars: individual 
(24 Sex, HD 200964) 
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1. INTRODUCTION 

The giant planets thus far discovered around other 
stars exhibit a wide variety of orbital characteristics that 
are very different from the properties of the planets in our 
Solar System. For example, exoplanets rarely reside in 
circular orbits and many of them have semimajor axes 1 
to 2 orders of magnitude smaller than those of Jupiter 
and Saturn. However, at least one characteristic of ex- 
oplanets appears to be shared with the constituents of 
the Solar System: planets often come in bunches. As 
Wright et al. 2009 recently showed, 28% of apparently 
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singleton exoplanets are later discovered to reside in sys- 
tems containing two or more components. As the pre- 
cision and time baselines of Doppler surveys increase, 
and as mor e planets are discovered from wide- field t ran- 
sit surveys (iBakos et al.ll2009bHCameron et alll2007D . di- 
rect imaging (iMarois et al.|l2 008: Ka las et alJl2008lj i. and 
microlensing (jGaudi et alJl2008b iGould et alJl2010f >. the 
currently measured multiplicity rate will likely prove to 
be a lower bound. 

Multiplanet systems are manifested in radial veloc- 
ity (RV) time series as either a single-planet Keplerian 
motion superimposed atop a partial longer -period orbit, 
or as multiple, time-resolved orbits (e.g. iWright et al.1 
2009). Those in the former category will gradually 
come into focus as the time baselines of Doppler sur- 
veys lengthen, and these "trend" systems are becoming 
recognized as promising direct-imaging targets. The lat- 
ter category, with their well-characterized orbits, are ex- 
tremely valuable for understanding the origins of plan- 
ets an d the evolution of their orbital architectures (|Fordl 
I2006bll . 

Observed deviations from pure Keplerian motions re- 
veal gravitational interactions among planets that serve 
as fossil records of past close encounters and/or conver- 
gent orbital migration (Malhotra 2002 IWu fc Goldreichl 
l2002f l. A prime exam ple is the system o f planets or- 
biting v Andromedae (|Butler et al.l Il999f) . iFord et all 
(2005) demonstrated through dynamical modeling that 
the current orbital configuration shows evidence of a vi- 
olent planet-planet scattering event in the distant past. 
Similarly, the two Jovian planets in the G1876 system 
currently reside in a mean-motion resonance (MMR) 



that may be a reflection of differential migration af- 
ter the planets ' form ation (jMarcv et al.ll200lllLeell2004l: 
lLaughlin et al.l [2005h . Thus, the discovery of exoplan- 
ets in MMRs is strong support for the inward orbital 
migration that is often invoked to explain the common 
existence of giant planets well within the canonical "snow 
line." 

Additionally, gravitational interactions among reso- 
nant planets can also place constraints on both the sys- 
tem inclination with respect to the sky, as well as mutual 
inclinations between the planets, and thereby remove the 
sin i ambiguity and provide absolut e measurements of th e 
planet masses (jRivera et al.ll2005t iCorreia et al.l I2010T) . 
Interactions observed in certain types of multiplanet sys- 
tems can reveal the interior structures of gas giant plan- 
ets in vivid detail. In the dramatic case of the system of 
planets around HAT-P-13, the inner planet transits its 
host star and experiences additional gravi tational per- 
turbations from an oute r planet near 1 AU (jBakos et al.l 
I2009al IWinn et al.ll2010D . Depending on the inclinations 
of the planets in the system, precise follow-up measure- 
ments may provide estimates of the tidal Love num- 
ber and Q value of the inner pla net to a higher pre- 
cision than is possible for Jupiter (|Batvgin et al.l 120091 : 
IMardlindl2010l) . 

We are conducting a Doppler survey of intermediate- 
mass subgiant stars at the Lick and Keck Observato- 
ries with the goal of understanding the influence stellar 
mass on the physical properties, orbital architectures and 
multiplicity rates of planetary systems. Our survey has 
resu lted in the discovery of 14 new singleton exoplan- 
ets (1 Johnson et al.|[2006U [20071 120081: [Peek et all l200l 



TABLE 1 

Stellar Parameters 



IBowler et al.ll2010t I Johnson et al.ll2010h . In this contri- 
bution we announce the discovery of two pairs of Jovian 
planets orbiting the subgiants 24 Sextanis (=HD 90043) 
and HD 200964. 

2. SPECTROSCOPIC OBSERVATIONS 

We began observations of 24 Sex and HD 200964 at 
Lick Observatory in 2004-2005 as part of our Doppler 
survey of intermediate-mass subgiants. Details of the 
survey, inclu ding target sele ction and observing s t rategy 
are given in [Jo hnson etUI (|2006al ). IPeek et all (|2009f ) 
and IBowler et alj (|2010t ). In 2007 we exp anded our sur- 
vey o f subgiants at Keck Observatory ([Johnson et al.1 
12010) and we added 24 Sex and HD 200964 to our Keck 
target list for additional monitoring. 

At Lick Observatory, the Shane 3 m and 0.6 m Coude 
Auxiliary Telescopes (CAT) feed the Hamilton spec- 
trometer (jVogtl 119871 ). and observations at Keck Ob- 
servatory were ob tained using the HIRES spectrometer 
(|Vogt et al.l I1994D . Doppler shifts are measured from 
ea ch o bservation using the iodin e cell method describe d 
bv iButler et~"aH (| 19961 ) (see also lMarcv fe ButleHli~992T) . 
A temperature-controlled Pyrex cell containing gaseous 
iodine is placed at the entrance slit of the spectrome- 
ter. The dense set of narrow molecular lines imprinted 
on each stellar spectrum from 5000 to 6000 A provides 
a robust wavelength scale for each observation, as well 
as information about the shape of the sp ectrometer's in- 
strumental response (jValenti et al.lll995l) . 

At Lick, typical exposure times of 60 minutes on the 
CAT and 5 minutes on the 3 m yield a signal-to-noise 
ratio (SNR) of wl20 at the center of the iodine region 



Parameter 


24 Sex a 


HD 200964 


V 


6.61 (0.04) 


6.64 (0.04) 


M v 


2.17 (0.06) 


2.35 (0.07) 


B-V 


0.92 (0.01) 


0.880 (0.009) 


Distance (pc) 


74.8 (4.9) 


68.4 (4.8) 


[Fe/H] 


-0.03 (0.04) 


-0.15 (0.04) 


Teff (K) 


5098 (44) 


5164 (44) 


Vrot sin i (km s ) 


2.77 (0.5) 


2.28 (0.5) 


log 9 


3.5 (0.1) 


3.6 (0.1) 


M, (Mo) 


1.54 (0.08) 


1.44 (0.09) 


R* (Rq) 


4.9 (0.08) 


4.3 (0.09) 


L* (Lq) 


14.6 (0.1) 


11.6 (0.4) 


Age (Gyr) 


2.7 (0.4) 


3.0 (0.6) 


logi?' HK 


-5.1 


-5.1 



a HD 90043 

(A = 5500 A), providing a velocity precision of 4.0- 
5.0 m s _1 . At Keck, typical spectra have SNR w 180 at 
5500 A, resulting in a velocity precision of 1.5-2.0 m s _1 . 
In addition to the internal, photon-limited uncertain- 
ties, the RV measurements include an additional noise 
term due to stellar "jitter" — velocity noise in excess 
of internal errors due to astrophysical sources such as 
pulsation and rotational modul ation of surface features 
ijSaaret all 1199ft I Wrightl I2005T ) . We therefore adopt a 
jitter val ue of 5 m s" 1 for our subgiants based on the es- 
timate of lJohnson et al.1 ([2010D . This jitter term is added 
in quadrature to the internal errors before determining 
the Keplerian orbital solutions. For the dynamical anal- 
ysis in §[5] we allow the jitter to vary as a free parameter 
in the fitting process. 

3. STELLAR PROPERTIES 

Atmospheric parameters of the target stars are esti- 
mated from iodine-free, "template" spectra using the 
LTE spectro scopic analysis package Spectroscopy Made 
Easy ( SME; Valcnti & Piskuno vl[1996D~ as descr i bed b y 
iValenti fe Fischerl ([2005D and iFischer fe Valentil (p005T ). 
To constrain the low surface gravit ies of the evolv e d star s 
we used the iterative scheme of IValenti et al.l (2009) , 
which ties the SME-derived value of log q to the gravity 
inferred from the Yonsei-Yale (Y 2 ; lYi et al.ll2004D stellar 
model grids. The analysis yields a best-fitting estimate 
of T e gf, logg, [Fe/H], and Vi-ot sin i. The properties of our 
targets from Lick and Keck are listed in the fourth edi- 
tion of the Spectroscopic Properties of Cool Stars Cat- 
alog (SPOCS IV.; Johnson et al. 2010, in prep). We 
adopt the SME parameter uncertainties described in the 
error analysis of Vale nti fc Fischerl (|2005f ) . 

The luminosity of each star is estimated from the ap- 
parent V-band magnitude and the parallax from Hip- 
parcos (|van Leeuwenl 120071). together wit h the b olomet- 
ric correction from IVa ndcnBcrg & Clem (|2003|) . From 
T e ff and luminosity, we determine the stellar mass, ra- 
dius, and an age estimate by associating those ob- 
served properties wit h a model fro m the Y 2 stellar 
interior calculations (|Yi et al.l 120041 ). We also mea- 
sure t he chromospheric e mission in the CallH&K line 
cores ([Wright et al.l 120041) . providing an S value on the 
Mt. Wilson system , which we convert to R' HK as per 
iNoves et al.l ([19841) . The stellar properties of the host 
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Fig. 1. — RV time series and two-Keplerian orbital solution for 
24 Sex. Top: Inner planet with the signal from the outer planet 
removed. The dashed line shows the best-fitting orbital solution 
for 24 Sexb. Middle: The outer planet with the inner planet re- 
moved. The dashed line shows the best-fitting orbital solution for 
24 Sexc. Bottom: The full RV time series, with the two-planet 
solution shown as a dashed line. 

stars are summarized in Table [TJ 

4. KEPLERIAN ORBITAL SOLUTIONS 

In this section we present the RV time series for both 
stars and the initial orbital analysis, which consists of the 
sum of two Keplerians without gravitational interaction. 
In § |6] we present the results of our Newtonian dynami- 
cal analysis for each system, which properly accounts for 
non-Keplerian motion due to gravitational interactions 
between the planets and host star. 

To search for the best-fitting, two-planet orbital solu- 
tion for each time seri es we use the partially-linearized 
technique presented by Wri ght fc Howard! (|2009f l. as im- 
plemented in the IDL software suite RVLIN. We es- 
timate the parameter uncertainties using a Markov- 
Chain Monte Carlo (MCMC) algorithm, as described by 
iBowleret all (|20Tol) . 

4.1. 24 Sextanis 

We obtained initial-epoch observations of 24 Sex at 
Lick Observatory in 2005 February, and since then we 
have obtained 50 RV measurements. After noticing time- 
correlated RV variations, we began additional monitoring 
at Keck Observatory in 2008 December, where we have 
obtained 24 additional RV measurements. The RVs from 
both observatories are listed in Table [H along with the 
Julian Dates (JD) of observation and the internal mea- 
surement uncertainties (without jitter). Fig. [T] shows the 
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Fig. 2. — RV time series and two-Keplerian orbital solution for 
HD 200964. Top: Inner planet with the signal from the outer planet 
removed. The dashed line shows the best-fitting orbital solution 
for HD 200964 b. Middle: The outer planet with the inner planet 
removed. The dashed line shows the best-fitting orbital solution for 
HD 200964 c. Bottom: The full RV time series, with the two-planet 
solution shown as a dashed line. 

RV time series from both observatories, where the error 
bars represent the quadrature sum of the internal errors 
a nd 5 m s -1 of jitter. 

iBowler et al.l (|2010D reported evidence of a two-planet 
system around 24 Sex, but the data at the time could 
not provide a unique orbital solution. An additional sea- 
son of observations has provided stronger constraints on 
the possible orbits of the two planets. We fitted a model 
consisting of two non- interacting planets and a 1.54 M© 
star orbiting their mutual center of mass. We find that 
two Keplerians provide an acceptable fit to the data with 
an rms scatter of 6.8 m s _1 and a reduced yXr/ = 1-14. 
The inner planet has a period of P = 455.2 ± 3.2 days, 
velocity semiamplitude K = 33.2 ± 1.6 m sA an d ec- 
centricity e = 0.184 ± 0.029. The outer planet has 
P = 910 ± 21 days, K = 23.5 ± 2.9 m sA, and e = 
0.412 ± 0.064. Together with our stellar mass estimate 
these spectroscopic orbital parameters yield semimajor 
axes {ajj,a c } = {1.41, 2.22} AU and minimum planet 
masses {Mb,M c } s'mi = {1.6,1.4} Mj up . A more de- 
tailed, dynamical analysis presented in § [6] revises this 
two-Keplerian solution under the constraint of long-term 
stability. 

4.2. HD 200964 

We began monitoring HD 200964 at Lick Observatory 
in 2007 July and have obtained 61 RV measurements. 
Time-correlated variations in the star's RVs prompted 



TABLE 2 
Orbital Parameters from Non-interacting Two-planet Solution 



Parameter 


24 Sexb b 


24 Sexc 


HD 200964 b 


HD 200964 b 


Period (d) 


455.2 (3.2) 


910 (21) 


630.6 (9.3) 


829 (21) 


T p b (JD) 


2454758 (30) 


2454941 (30) 


2454916 (30) 


2455029 (130) 


Eccentricity 


0.184 (0.029) 


0.412 (0.064) 


0.111 (0.030) 


0.113 (0.076) 


K (m s" 1 ) 


33.2 (1.6) 


23.5 (2.9) 


35.2 (2.7) 


22.1 (2.3) 


u> (deg) 


227 (20) 


172.0 (9) 


223 (20) 


301 (50) 


M P sin i (M Jup ) 


1.6 (0.2) 


1.4 (0.2) 


1.9 (0.2) 


1.3 (0.2) 


a (AU) 


1.41 (0.03) 


2 


22 (0.06) 


1.71 (0.04) 


2.03 (0.06) 


Lick rms (m s _1 ) 


7.7 






7.6 




Keck rms (m s —1 ) 


1.8 






5.3 




Jitter (m s _1 ) 


5.0 






5.0 




^2 


1.14 






1.15 




N oba Lick 


50 






61 




N obs Keck 


21 






35 





a HD 90043 

b Time of periastron passage. 

additional monitoring at Keck Observatory where we 
have obtained 35 measurements since 2007 October. The 
RVs from both observatories are listed in Table [SJ along 
with the Julian Dates (JD) of observation and the inter- 
nal measurement uncertainties (without jitter). Fig. [2] 
shows the RV time series from both observatories, where 
the error bars represent the quadrature sum of the inter- 
nal errors and 5 m s _1 of j itter. 

As is the case for 24 Sex. lBowler et al. (2010) reported 
evidence of a two-planet system around HD 200964, 
and an additional season of observations has provided 
stronger constraints on the possible orbits of the plan- 
ets in the system. We find that a two-Keplerian model 
provides an acceptable fit to the data with an rms scat- 
ter of 6.8 m s _1 and a reduced \/xf — 1-14. The in- 
ner planet has a period of P = 630.6 ± 9.3 days, ve- 
locity semiamplitude K = 35.2 ± 2.7 m s -1 , and ec- 
centricity e = 0.111 ± 0.030. The outer planet has 
P = 829 ± 21 days, K = 22.1 ± 2.3 m s" 1 , and e = 
0.113 ± 0.076. Together with our stellar mass estimate 
these spectroscopic orbital parameters yield semimajor 
axes {af,,a c } — {1.71, 2.03} AU and minimum planet 
masses {M},,M C } sini = {1.9,1.3} Mj up . A more in- 
depth dynamical analysis presented in § |6] revises this 
two-Keplerian solution under the constraint of long-term 
stability. 

5. PHOTOMETRIC MONITORING 

In addition to the radial velocities presented in § U we 
also acquired photometric measurements of both 24 Sex 
and HD 200964 with the T12 0.80 m automatic photo- 
metric telescope (APT) at Fairborn Observatory. The 
T12 APT and its two-channel photometer measure count 
rates simultaneously through Stromgren b and y filters. 
T12 is ess entially identi cal to the T8 0.80 m APT de- 
scribed in iHenrvi (|1999[ ). Each program star (P) was 
observed differentially with respect to two nearby com- 
parison stars (CI and C2) (see Table [6]). The differen- 
tial magnitudes P - CI, P - C2, and C2 - CI were 
computed from each set of differential measures. All dif- 
ferential magnitudes with internal standard deviations 
greater than 0.01 mag were rejected to eliminate observa- 
tions taken under non-photometric conditions. The sur- 
viving observations were corrected for differential extinc- 



tion with nightly extinction coefficients and transformed 
to the Stromgren photometric system with yearly-mean 
transformation coefficients. We averaged the b and y 
measurements of each star into a single (b + y)/2 "pass- 
band" (which we designate by in Table [6]) to improve the 
precision of the brightness measurements. Typical pre- 
cision of a single (b + y)/2 observation, as measured for 
pairs of constant stars, is ~0. 0015-0. 0025 mag on good 
p hotometric nights. 

lOuelozet alT(|2001[ ) and iPaulson et all (|2004D have 
demonstrated how rotational modulation in the visibility 
of star spots on active stars can result in periodic radial 
velocity variations that mimic the presence of a planetary 
companion. Thus, the precise APT brightness measure- 
ments are valuable for distinguishing between activity- 
related radial velocity changes and true reflex motion of 
a star caused by a planet. 

Photometric results for 24 Sex and HD 200964 are 
given in Table [6] Columns 7-9 give the standard devia- 
tions of the P-Cl, P-C2, and C2-C1 differential mag- 
nitudes in the (b + y)/2 passbands. All of the standard 
deviations are small and within the range of measure- 
ment precision with the T12 APT. We also performed 
periodogram analyses on all the data sets and found no 
significant periodicity between 0.03 and 100 days that 
might be the signature of stellar rotation. Our data sets 
are not sufficiently long to test for variability on the four 
planetary periods from Table [2] but we expect any such 
variability to be very small. The photometric stability of 
both 24 Sex and HD 200964 and the long-term coherency 
of the observed radial velocity variations provide strong 
support for the existence of all four new planets. 

6. DYNAMICAL INTERACTIONS 

Our best-fitting double-Keplerian fits have two 
Jupiter-mass planets with periods near the 2:1 commen- 
surability for 24 Sex and near 4:3 for HD 200964. The 
proximity of the two planets implies strong gravitational 
interactions, which limits the number of possible orbits 
to those that allow the two planets to remain stable over 
the lifetime of the star. To test the long-term stability 
of the various orbital configurations that are consistent 
with the data we performed a series of numerical inte- 
grations as described in the following sections. 



6.1. Methodology for MCMC Analysis incorporating 
N-Body Integrations 

We analyze the radial veloci ty ob s ervatio ns us ing a 
Bayesia n framework following iFordl (|2005f) and iFordl 
l|2006af) . We assume priors that are uniform in log- 
arthimic intervals of orbital period, eccentricity, argu- 
ment of pericenter, mean anomaly at epoch, and the ve- 
locity zero-point. For the velocity amplitude (K) and 
stellar jitter (<Jj), we adopt a prior of the form p(x) — 
(x + x )~ 1 [log(l + x/x )]~~ 1 , with K a = g 1p = lm/s . 
For a discussion of priors, see iFord fe Gregory! (|2007j ). 
The likelihood for radial velocity terms assumes that each 
radial velocity observation (vi) is independent and nor- 
mally distributed about the true radial velocity with a 
variance of of + cr? where o~i is the published measure- 
ment uncertainty and o~j is a jitter parameter that ac- 
counts for additional scatter due to stellar variability, in- 
strumental errors and/or inaccuracies in the model (i.e., 
neglecting planet-planet interactions or additional, low- 
amplitude planet signals). 

In our initial phase of analysis, we use an MCMC 
method based upon Keplerian orbit fitt ing to calcu late a 
sample from the posterior distribution ([Ford 2006a| . We 
calculate multiple Markov chains, each with ~ 2 x 10 8 
states. We discard the first half of the chains and calcu- 
late Gelman-Rubin test statistics for each model param- 
eter and several ancillary variables. We find no indica- 
tions of non-convergence. Thus, we randomly choose a 
subsample (^16, 000 samples) from the posterior distri- 
bution for further investigation. 

Next, we use this subsample as the basis for a much 
more computationally demanding analysis that uses fully 
self-consistent n-body integrations to account for planet- 
planet interactions when modeling the RV observations. 
We again perform a Bayesian analysis, but replace the 
standard MCMC algorithm with a Differential Evolu- 
tion Markov chain Monte Carlo (DEMCM C) algorithm 
(jter Braakll2006HVeras fc Fo7l2"009l 120101 ) . In the DEM- 
CMC algorithm each state of the Markov chain is an 
ensemble of orbital solutions. The candidate transition 
probability function is based on the orbital parameters in 
the current ensemble, allowing the DEMCMC algorithm 
to sample more efficiently from high-dimensional param- 
eter spaces that have strong correlations between model 
parameters. More details of this DEMCMC algorithm 
and associated tests of its accuracy will be presented in 
a forthcoming paper, (Nelson et al., 2011, in prep.) 

The priors for the model parameters are the same as 
those of the MCMC simulations. The initial conditions 
of the n-body simulations are calculated by converting 
between Keplerian and Cartesian coordinates. In this 
paper, we consider only coplanar two-planet systems. 

For the n-body integrations, we use a time symmet- 
ric 4th order Hermite integrator that has been optimized 
for planetary systems (Kokubo et al. 1998). We extract 
the radial velocity of the star (in the barycentric frame) 
at each of the observation times for comparison to RV 
data. During the DEMCMC analysis, we also impose 
the constraint of short-term (100 years) orbital stability. 
We check whether the planetary semimajor axes remain 
within a factor of 50% of their starting value, and that no 
close-approaches occur within 0.1 x the semimajor axis 
during the 100 year n-body integration. Any systems 



failing these tests are rejected as unstable (regardless of 
the quality of the fit to RV data). Thus, the DEMCMC 
simulations avoid orbital solutions that are violently un- 
stable. In our DEMCMC simulations, this process is 
repeated for 10,000 generations, each of which contains 
16,000 systems, for a total of ~ 10 8 n-body integrations 
in each DEMCMC simulation. 

Since the DEMCMC simulations only require stability 
for ~ 100 years, the orbital solutions in the final gen- 
eration may or may not be stable for longer timescales. 
Since nearly all of these systems are strongly interact- 
ing, we take this final generation (16k systems) and de- 
mand that they also be stable (according to the same 
criteria above) over the course of a 10 7 year integration, 
performed using the hybrid Bulirsch-Stoer / Symplectic 
integrator Mercury (Chambers 1999). Only the orbital 
solutions which are stable over the course of this long- 
term n-body integration are regarded as being accept- 
able solutions. More details on the dynamical analysis 
performed and the results obtained will be presented in 
a companion paper, (Payne et al, 2010, in prep.). 

6.2. Numerical Integrations for 24 Sex 

We find that the n-body DEMCMC routine results for 
24 Sex are concentrated around solutions with {Pb, P c } ~ 
{450,900} days, i.e. close to the 2:1 resonance, and 
that they span a significantly smaller range of parame- 
ter space than do the double-Keplerian fits. The n-body 
RV fitting aspect of the DEMCMC routine acts to shrink 
the parameter space, whereas the stability requirement 
in the routine has the effect of shifting the DEMCMC 
solutions towards the 2:1 period ratio mark, primarily 
by favoring lower periods for the inner planet. We illus- 
trate in Fig. [3]the difference in the end-points of the two 
analyses, showing a scatter plot of the planetary periods 
at the end of both the Keplerian analysis and the n-body 
analysis. 

We have applied the DEMCMC method described in 
§ 16.11 multiple times using different values for the initial 
ensemble of orbital solutions. Each time, these simula- 
tions reached qualitatively similar results to those of Fig. 
[31 but quantitatively there are signs that the method has 
not fully converged. Therefore, we do not interpret the 
results as a precise estimate of the posterior probabil- 
ity distribution. Instead, we use these results to demon- 
strate that there exist orbital solutions that are both sta- 
ble and consistent with the Doppler observations. Most 
importantly, all our simulations indicate that the poste- 
rior distribution is concentrated at solutions with a ratio 
of orbital periods between ~ 1.9 and ^2.1 and very close 
to the 2:1 MMR. Thus, we conclude that the current ob- 
servations strongly favor orbital solutions in (or at least 
near) the 2:1 MMR. 

Next, we performed long-term stability tests for the 
16,000 orbital solutions in the final generation of the n- 
body DEMCMC analysis. We find that ~ 90% of the 
orbital solutions are unstable over the course of 10 7 years 
of integration. We consider the remaining ~ 1 , 000 (of 
the 16,000 orbital solutions) that are stable for 10 7 years 
to be plausible orbital solutions given both the RV data 
for 24 Sex and the requirement of long-term dynamical 
stability. 

We show in Fig. [3] that the stable systems remain 



TABLE 3 
Orbital Parameters from Stable, Interacting N-Body Two-planet Solution 



Parameter 


24 Sexb b 


24 Sex c 


HD 200964 b 


HD 200964 c 


Period (d) 


452.8 (tH) 


883.0 (til) 


613.8 (tii) 


825.0 (til) 


T p b (JD) 


2454762 (tf 7 l 3 ) 


2454930 ( + gg 9 5 9 ) 


2454900 ( + g| 5 3 5 ) 


2455000 (t!il) 


Eccentricity 


o.o9 nit) 


0-29 (t° .H) 


0.04 (+°;°*) 


0-181 (±8:23) 


K (m s" 1 ) 


40.0 (t 4 7 9 s ) 


14-5 (tH) 


34.5 (i? : J) 


15-42 (till) 


w (deg) 


9-2 (t^i) 


220.5 (_32o.g) 


288.0 (±fli° 9 ) 


182.6 (±^;I) 


M P sin i (A/jup) 


i-99 (±g;ll) 


0.86 (t°;i) 


1-85 (lo'-.os) 


0-90 (±°;») 


a (AU) 


1-333 (t° .Z) 


2.08 (+°;«) 


1-601 (°J 2 02 ) 


1 qz. /+0.008x 
J-WO i-0.005^ 


Jitter (m s — ) 


9.9 (i?: 9 2 ) 




8-23 (±°;11) 





strongly clustered around the 2:1 period commensura- 
bility region. Taking all of these long-term stable sys- 
tems into account, we find that the data indicate that 
the inner planet has a period P& = 452.8 t.4'5 days, 
semimajor axis a& = 1.333to'oo9 AU, eccentricity e b = 

h0.26 



0.09lQQg and mass M b sin i 

outer planet has a period P c = 883.0^4 days, semimajor 



1.9912;^ M Jup , while the 



axis a r 



+0 ° 5 2 AU, eccentricity e c = 0.29±g;J| and 
0.86t°;22 M Ju P - We detail all the fitted 



- = 2.08 
mass Af c sin i = 0.86^0:22 M, 
parameters from this analysis in Table [3) 

When we examine specific systems in detail, we find 
that the pericenter of the outer planet overlaps the loca- 
tion of the pericenter of the inner planet, meaning that 
the planets therefore experience detectable gravitational 
interactions. The systems typically undergo large eccen- 
tricity oscillations (0.24 < e b < 0.47, 0.5 < e c < 0.6) 
over a period of ~ 800 years, while the semimajor axis 
varies on a shorter (~ 25 year) timescale, with a range 
1.29 <a b < 1.31 and 2.04 < a c < 2.10. 

We note that the PDFs from the DEMCMC analysis 
are bimodal. To investigate the cause of the biomodal- 
ity, we performed a similar n-body+DEMCMC analysis 
without the requirement of short-term stability. The re- 
sulting PDFs have a single broader mode consistent with 
the output of Keplerian MCMC analysis. We conclude 
that the biomodal nature of the PDF from the dynamical 
analysis is most likely the result of demanding short-term 
stability. 

The overlap of the pericenters suggests that a MMR 
is needed to stabilize the system over long timescales by 
preventing close encounters. However, an analysis of the 
resonant angles #2,1 = 2A C — lAj, — liz'2,1 (where A is 
the mean longitude and w is the longitude of pericenter) 
suggests that, for all the long-term-stable systems, the 
planets circulate, rather than librate, i.e. all of the sys- 
tems are observed to have angular ranges for #2.1 ~ 360°. 
We therefore caution that the current state of the obser- 
vations and dynamical analysis cannot confirm whether 
the system is truly in a resonant state, and hence fur- 
ther work will be required. A more detailed investiga- 
tion (Payne et al. 2010 in prep) will probe the nature of 
these dynamical interactions in 24 Sex and HD 200964 in 
greater detail. 

Finally, we note that as discussed in § 16.11 jitter is al- 
lowed to vary during the DEMCMC procedure. As such 
we find a best fit value for the jitter in the same manner 



as we do for the various planetary orbital parameters. 
From the 24 Sex analysis, we find a large jitter value of 
9>9ii 2 m s_1 - This is substantially larger than the value 
of 5 m s" 1 assumed in § |4j indicating that 24 Sex has 
intrinsic RV variability at the high end of the observed 
distribution, or that there are other unmodeled sources 
of variability such as an additional planet in the system. 
However, a periodogram of the residuals shows no sig- 
nificant power abov e the noise. Using the approach of 
I Bowler et al.l (|2010D . we can rule out additional compan- 
ions with masses greater than 0.3 Mj up out to 1 AU and 
"hot" planets with masses Mp sin i < 0.1 Mj up within 
0.1 AU. Additional monitoring at Keck and a more in- 
depth dynamical analysis will help clarify the situation. 

6.3. Numerical Integrations for HD 200964 

The results of our n-body DEMCMC simulations for 
HD 200964 are again more tightly confined than the re- 
sults of the double-Keplerian analysis. As an example 
of this we illustrate in Fig. |4] the difference in the end- 
points of the two analyses, showing a scatter plot of the 
planetary periods at the end of both the Keplerian analy- 
sis and the n-body analysis, as well as histograms for the 
same data. Again, we note that the shift to an N-Body 
analysis acts to concentrate the preferred region into a 
smaller area, while the stability requirements act to shift 
the favored parameters closer toward the 4:3 period com- 
mensurability 

As for the case of 24 Sex, we have applied our DEM- 
CMC method multiple times using numerous randomized 
initial conditions. Again, we find qualitatively similar re- 
sults from all sets of the DEMCMC analysis, but with 
some indication that the results have not truly converged. 
As such, we interpret these results as demonstrative that 
stable systems exist which can fit the RV data, and more- 
over, the extremely narrow range of period ratios favored 
by the analysis (~ 1.32 to ~ 1.36) shows that the current 
observations strongly favor orbital solutions in or close to 
the 4:3 MMR. 

Next, we test the long-term orbital stability of the or- 
bital solutions identified by the DEMCMC analysis. For 
HD 200964, > 90% of the systems are clearly unstable 
during a 10 7 year integration (i.e., they experience a col- 
lision or a change in semimajor axes of more than 50%.) 
We find that ~ 1,000 of the 16,000 simulations remain 
stable for 10 7 years. We interpret these as plausible or- 
bital solutions consistent with the RV observations of HD 
200964. 
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Fig. 3. — Top: Comparison of the end-points of the double- 
Keplerian MCMC analysis (shaded contours), and the n-body 
DEMCMC analysis with the long-term stable constraint (black, 
solid contours) for 24 Sex. The contours show the 68.2%, 80% and 
95% confidence intervals (inner to outer). To guide the eye, we 
overlay diagonal lines to indicate 1.75:1, 2:1 & 2.2:1 period ratios. 
The MCMC posterior sample covers a large range of parameter 
space and is consistent with a 2:1 period ratio. In contrast, the 
DEMCMC results are limited to a smaller area of parameter space 
lying on, or close to, the 2:1 period ratio. Bottom: Histograms for 
the same planetary period data from 24 Sex, comparing the results 
from the Keplerian MCMC analysis (black), N-Body DEMCMC 
analysis (dark gray). When compared to the basic MCMC results, 
the DEMCMC results for the inner planet are shifted more towards 
lower periods, while the outer planet shifts towards higher periods, 
making the overall period ratio converge more closely towards 2:1. 

As illustrated in Fig. [H the stable systems occupy 
a region of parameter space corresponding to a region 
of paramter space near the 4:3 period ratio. Taking 
only these long-term stable systems into account, we find 
that the inner planet has a period Pb = 613.8^'| days, 
semimajor axis a& = 1.601°_q 2 02 AU, eccentricity e b = 

0.04±o.u2 and mass M b sini = lMto.os M Ju P , while the 
outer planet has a period P c — 825. 0J^ '\ days, semima- 

+0.024 
-0.017 



j or axis a c 



i nc+0.008 
1 - yd -0.005 



AU, eccentricity e c = 0.181 



and mass M c sini = 0.90Ioog M Jup . We detail all the 
fitted parameters from this analysis in Table [3] 

We find that most of the stable planetary orbits over- 
lap, producing strongly interacting systems, resulting in 
significant oscillations in the the semimajor axes and 
the eccentricities of both planets. An example sta- 
ble solution exhibits observable eccentricity oscillations 
(0.03 < e b < 0.1 and 0.13 < e c < 0.18) on an ap- 
proximately 250-year timescale and semimajor axis os- 
cillations (1.57 < a,b < 1.58 and 1.9 < a c < 1.93) on 
an approximately 60 year timescale. Again, as in the 
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Fig. 4. — Top: Comparison of the end-points of the double- 
Kcplcrian MCMC analysis (shaded contours), and the n-body 
DEMCMC analysis with the long-term stable constraint (black, 
solid contours) for 24 Sex. The contours show the 68.2%, 80% and 
95% confidence intervals (inner to outer). To guide the eye, we 
overlay diagonal lines to indicate 1.15:1, 4:3 & 1.45:1 period ratios. 
The MCMC posterior sample covers a large range of parameter 
space and is consistent with a 4:3 period ratio. In contrast, the 
DEMCMC results are limited to a smaller area of parameter space 
straddling the 4:3 period ratio. Bottom: Histograms for the two 
planetary periods in the HD 200964 system, comparing the results 
from the Keplerian MCMC analysis (black) and N-Body DEM- 
CMC analysis (dark gray). When compared to the basic MCMC 
results, the DEMCMC results in a decrease in the fitted period 
of the inner planet, making the overall period ratio converge more 
closely towards 4:3. 

24 Sex system, all of the long-term stable systems ex- 
amined here for HD 200964 seem to be circulating rather 
librating. As in the 24 Sex analysis, we find a jitter value 
of 8.23lQ'ggm s _1 , which is larger than the empirical jit- 
ter estimate of Uohnson et al.l ((2010). 



7. SUMMARY AND DISCUSSION 

Our RV measurements of the intermediate-mass sub- 
giants 24 Sex and HD 200964 (A/* = 1.54 and 1.44 M©, 
respectively) reveal the presence of a pair of giant plan- 
ets around each star. Our orbital analysis indicates 
that most, if not all, of the dynamically stable solutions 
include crossing orbits, suggesting that each system is 
locked in a mean motion resonance that prevents close 
encounters and provides long-term stability. The planets 
in the 24 Sex system likely have a period ratio near 2:1, 
while the HD 200964 system is even more tightly packed 
with a period ratio close to 4:3. 

In both the 24 Sex and HD 200964 systems, the plan- 
ets reside well within the so-called snow line, beyond 
which volatiles in the protoplanetary disk can condense 



to provide the raw materials for protoplanetary core 
growth. For a pre- main-sequence, 1.5 M© star the snow 
line is lo cated beyond 2-3 AU acco rding to the esti- 
mates of Kennedy fe Kenvonl (|2008l ) for realistic disk 
models including irradiation and accretional heating. 
It is therefore likely that the planets around 24 Sex 
and HD 200964 formed at larger semimajor axes and 
subsequently experienced inw ard orbital migration (See 
iPapaloizou fc Terq ucm 2006, for a review of migration 
theory). 

Unless the planets around 24 Sex and HD 200964 
formed with period ratios close to their current values, 
the two planets in each system likely migrated through 
disk interactions at convergent rates until they became 
trapped in an MMR, with the strong 2:1 resonance be- 
ing t he most likely endpoint of such differe ntial migra- 
tion (jKlevll200Clt iNelson fc Papaloizoull2002h . The rarity 
of planets discovered with period ratios smalle r than 2:1 
accord s well with the dynamical simulations of lLee et al.l 
(2009). In their simulation they considered the formation 
of two giant planets in a protoplanetary disk with initial 
period ratios just outside of 2:1 and final orbital config- 
urations determined by the initial conditions and details 
of the planet-planet-disk interactions. From their ensem- 
ble of simulated planetary systems they found that only 
3% attain period ratios closer than 2:1, and none ended 
in stable configurations closer than 3:2. 

The ability of planets to reach a resonance tighter 
than 3:2 is dependent upon a number of factors includ- 
ing the initial separation of the planets, disk viscosity, 
planet masses, remaining disk mass and size of the gaps 
open e d by the planets (Malhotra 1993; Haghigl npourl 
[19991: iBrTden et al.l 120001: ISnelkrove et al.l 12001ft . One 
key variable is the convergent migration rate, which if 
fast enough can carry the planets past the deep 2:1 MMR 
into tighter commcnsurabilities. The 24 Sex system is 
near the 2:1 resonance, which may be a reflection of a 
slower migration process compared to that which led to 
the extremely tight 4:3 configuration in the otherwise 
v ery similar HD 200964 system . 

Papaloiz ou fc Szuszkiewica ()2010h explored rapid mi- 
gration scenarios leading to the attainment of high-order 



MMRs by low-mass planets migrating withing a gas disk. 
For the terrestrial planets they considered, MMRs with 
degrees as high as 8:7 and 11:10 were achieved for mi- 
gration timescales of order 10 3 years. However, to test 
whether such conditions can lead to high-order stable 
MMRs for the planets in the HD 200964 system, and 
the 2:1 MMR seen in the 24 Sex system, hydrodynamic 
considerations ne ed to be incorporated into the models. 
iRein et ail (|2010D investigated the formation and evolu- 
tion of the gas giants orbiting HD 45364 and found plausi- 
ble models for the attainment of the 3:2 MMR observed 
in that system. Simulations of this nature are beyond 
the scope of the present work and will be presented in 
a future contribution (Payne et al. 2010, in prep). For 
now, it is clear that just like the resonant planetary sys- 
tems discovered before them, the 24 Sex and HD 200964 
systems pose interesting challenges to theories of planet 
formation and orbital evolution. 
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TABLE 4 
Radial Velocities for 24 Sextanis 



JD 


RV 


Uncertainty 


Telescope 


-2440000 


(m s _1 ) 


(m s~ 1 ) 




13405.936 


-34.98 


4.24 


L 


13435.836 


-37.40 


4.05 


L 


13502.774 


3.04 


3.76 


L 


13517.768 


8.07 


3.90 


L 


13719.011 


-13.68 


3.46 


L 


14131.051 


38.40 


4.17 


L 


14133.931 


37.92 


3.49 


L 


14146.787 


28.17 


4.15 


L 


14150.947 


34.70 


6.03 


L 


14168.823 


28.25 


4.33 


L 


14169.860 


30.86 


3.65 


L 


14178.699 


36.12 


4.10 


L 


14216.839 


-6.66 


4.22 


L 


14218.709 


-8.74 


6.44 


L 


14232.713 


-4.29 


4.05 


L 


14254.688 


-24.93 


6.01 


L 


14255.689 


-33.28 


4.29 


L 


14266.699 


-25.63 


4.77 


L 


14269.695 


-28.99 


5.87 


L 


14274.719 


-31.09 


4.86 


L 


14378.045 


-24.52 


6.79 


L 


14402.054 


2.53 


5.26 


L 


14433.046 


24.41 


5.29 


L 


14445.938 


4.46 


3.08 


L 


14458.050 


4.53 


4.31 


L 


14482.959 


11.65 


3.89 


L 


14504.855 


8.92 


3.36 


L 


14525.946 


25.65 


4.37 


L 


14548.900 


22.38 


3.29 


L 


14572.872 


13.31 


4.54 


L 


14573.724 


16.22 


3.77 


L 


14617.764 


11.34 


5.29 


L 


14620.705 


-1.69 


3.78 


L 


14621.716 


11.26 


5.22 


L 


14627.708 


11.39 


3.87 


L 


14650.692 


-2.99 


6.79 


L 


14762.024 


-24.44 


7.00 


L 


14786.078 


-13.05 


6.03 


L 


14806.150 


-22.03 


1.98 


K 


14807.075 


-18.78 


1.83 


K 


14808.057 


-13.83 


1.96 


K 


14809.069 


-21.81 


1.88 


K 


14810.149 


-15.43 


1.88 


K 


14811.122 


-12.80 


1.65 


K 


14847.038 


3.99 


1.80 


K 


14847.048 


18.53 


3.83 


L 


14904.830 


52.37 


4.91 


L 


14929.814 


59.40 


1.63 


K 


14955.878 


69.74 


1.60 


K 


14963.915 


68.12 


1.49 


K 


14978.706 


55.29 


4.31 


L 


14983.770 


64.81 


1.59 


K 


14984.827 


48.72 


1.53 


K 


14987.759 


63.15 


1.54 


K 


14988.750 


63.85 


1.60 


K 


15002.711 


55.03 


6.40 


L 


15014.741 


57.09 


1.72 


K 


15015.750 


61.18 


1.58 


K 


15151.052 


-23.66 


6.04 


L 


15164.112 


-21.86 


1.83 


K 


15172.150 


-26.29 


2.05 


K 


15173.097 


-31.15 


1.93 


K 


15189.119 


-40.47 


1.85 


K 


15229.007 


-30.27 


3.76 


L 


15241.871 


-40.78 


4.23 


L 


15242.872 


-36.26 


3.75 


L 


15243.855 


-43.99 


5.27 


L 


15260.941 


-26.13 


1.65 


K 


15271.797 


-27.09 


3.77 


L 


15272.789 


-11.58 


6.00 


L 


15273.750 


-25.81 


3.66 


L 


15285.854 


-12.17 


1.75 


K 


15311.810 


-5.48 


1.43 


K 


15312.791 


2.81 


1.56 


K 
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TABLE 5 
Radial Velocities for HD 200964 



JD 


RV 


Uncertainty 


Telescope 


-2440000 


(m s _1 ) 


(m s~ 1 ) 




13213.895 


11.40 


4.20 


L 


13255.775 


29.55 


3.72 


L 


13327.604 


36.56 


3.70 


L 


13576.944 


-20.58 


3.98 


L 


13619.810 


-36.28 


3.71 


L 


13641.795 


-27.30 


3.61 


L 


13669.629 


-18.07 


3.55 


L 


13710.605 


-4.99 


3.77 


L 


13894.977 


-3.08 


3.71 


L 


13895.921 


6.24 


3.30 


L 


13896.963 


5.08 


3.47 


L 


13921.953 


9.72 


3.59 


L 


13959.788 


7.02 


3.96 


L 


14216.976 


-19.74 


4.09 


L 


14232.936 


-26.06 


4.83 


L 


14244.907 


-27.39 


5.19 


L 


14254.960 


-16.97 


3.70 


L 


14274.916 


-7.88 


3.72 


L 


14288.888 


-3.15 


3.49 


L 


14304.877 


8.18 


4.89 


L 


14309.845 


-1.57 


4.28 


L 


14336.881 


16.00 


4.49 


L 


14377.773 


10.58 


5.74 


L 


14399.752 


33.36 


1.62 


K 


14401.749 


41.83 


5.67 


L 


14405.700 


39.39 


4.02 


L 


14427.588 


35.30 


4.29 


L 


14427.757 


25.14 


1.18 


K 


14429.621 


49.55 


4.41 


L 


14432.653 


30.85 


4.89 


L 


14445.624 


23.13 


4.22 


L 


14536.053 


20.85 


5.26 


L 


14551.019 


25.15 


7.90 


L 


14585.990 


22.22 


4.51 


L 


14603.126 


14.71 


1.39 


K 


14604.012 


18.01 


1.58 


K 


14612.988 


-7.49 


6.44 


L 


14621.965 


3.21 


5.14 


L 


14622.995 


-5.37 


4.16 


L 


14634.079 


10.67 


1.52 


K 


14640.923 


-1.39 


3.34 


L 


14650.960 


-13.59 


3.79 


L 


14656.902 


0.00 


7.81 


L 


14674.916 


-16.16 


1.39 


K 


14675.850 


-16.06 


4.87 


L 


14676.879 


-25.35 


4.49 


L 


14677.922 


-28.67 


3.69 


L 


14683.850 


-22.03 


4.24 


L 


14699.804 


-14.57 


5.28 


L 


14734.713 


-29.88 


3.82 


L 


14737.790 


-29.53 


5.34 


L 


14738.751 


-36.28 


4.54 


L 


14747.711 


-59.70 


3.80 


L 


14766.659 


-48.87 


3.96 


L 


14778.803 


-45.21 


1.61 


K 


14785.605 


-51.33 


4.06 


L 


14790.734 


-47.17 


1.34 


K 


14807.789 


-46.94 


1.60 


K 


14935.136 


-18.39 


1.16 


K 


14941.012 


-18.06 


6.19 


L 


14941.989 


-24.98 


4.38 


L 


14956.097 


-14.88 


1.50 


K 


14964.120 


-5.85 


1.33 


K 


14978.935 


-3.81 


5.48 


L 


14984.070 


6.28 


1.43 


K 


14985.095 


4.20 


1.62 


K 


14986.112 


-1.27 


1.71 


K 


14987.129 


3.06 


1.39 


K 


14989.069 


6.64 


1.33 


K 


15002.947 


32.46 


5.34 


L 


15014.972 


22.36 


1.54 


K 


15015.957 


11.56 


1.44 


K 


15027.007 


22.52 


1.55 


K 


15042.973 


33.72 


1.69 


K 


15049.001 


42.16 


1.62 


K 


15060.816 


45.09 


4.29 


L 


15061.819 


45.31 


3.69 


L 


15062.829 


47.93 


3.69 


L 


15075.078 


45.22 


1.68 


K 


15076.067 


43.59 


1.75 


K 
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TABLE 6 
Summary of Photometric Observations From Fairborn Observatory 



Program 
Star 
(1) 



Comparison 
Star 1 

(2) 



Comparison 
Star 2 

(3) 



Date Range 
(HJD - 2,400,000) 

(4) 



Duration 

(days) 

(5) 



N obs 
(6) 



"(P-Cl) by 

(mag) 

(7) 



°{P - C2) by 

(mag) 

(8) 



a(C2 - Cl) by 

(mag) 

(9) 



Variability 
(10) 



24 Sex 
HD 200964 



HD 87974 
HD 201507 



HD 89734 
HD 201982 



54438-55333 
54578-55167 



895 



243 
100 



0.0022 
0.0015 



0.0020 
0.0014 



0.0020 
0.0015 



Constant 
Constant 



